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ABSTRACT 

During the first 20 seconds of its life, the enormous neutrino luminosity of a 
neutron star drives appreciable mass loss from its surface. This neutrino-driven 
wind has been previously identified as a likely site for the r-process. Qian & 
Woosley (1996) have derived, both analytically and numerically, the physical 
conditions relevant for heavy element synthesis in the wind. These conditions 
include the entropy (S), the electron fraction (Y e ), the dynamic time scale, 
and the mass loss rate. Here we explore the implications of these conditions 
for nucleosynthesis. We find that the standard wind models derived in that 
paper are inadequate to make the r-process, though they do produce some rare 
species above the iron group. We further determine the general restrictions on 
the entropy, the electron fraction, and the dynamic time scale that are required 
to make the r-process. In particular, we derive from nuclear reaction network 
calculations the conditions required to give a sufficient neutron-to-seed ratio 
for production of the platinum peak. These conditions range from Y e ps 0.2 
and S < 100 per baryon for reasonable dynamic time scales of ~ 0.001-0.1 
s, to Y e ks 0.4-0.495 and S > 400 per baryon for a dynamic time scale of 
~ 0.1 s. These conditions are also derived analytically to illustrate the physics 
determining the neutron-to-seed ratio. 

Subject headings: elementary particles — nuclear reactions, nucleosynthesis, 
abundances — supernovae: general 
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1. Introduction 

The origin in nature of the neutron-rich isotopes heavier than the iron group remains 
uncertain, but rapid progress in both theory and observation has been made in recent years. 
For example, analysis of the metal deficient star CS 22892-052 ([Fe/H] ~ —3.1) by Sneden 
et al. (1996) shows striking evidence that the solar abundance pattern for the r-process 
isotopes exists even at very early times in our Galaxy, and suggests an origin for these 
nuclei quite distinct from the s-process. The fact that the r-process abundance pattern in 
this very metal deficient star is so strikingly solar across the entire range 56 < Z < 76 also 
suggests that this pattern is generic and reflects the conditions constantly obtained in a 
unique astrophysical environment responsible for the r-process. However, we note that the 
current observational data for Z < 56 are not so conclusive. 



On the theoretical front, many recent calculations (e.g., [Woosley fc Hoffman 1992 



Howard et al. 1993| ; [Takahashi, Witti, fc Janka 1994| ; [Woosley et al. 1994Q have shown 



that a promising site for r-process nucleosynthesis is the neutrino-driven wind blowing 
from a nascent neutron star following a core-collapse-driven supernova. This r-process site 
has some attractive features. First of all, the r-process would be primary, in accordance 
with the observation of Sneden et al. 19961 (see also Cowan et al. 1996 ; Mathews fc 



Uowan 1990 ; |Cowan, Thielemann, fc Truran 199 1| ). Furthermore, the total mass loss in the 



wind approximately accounts for the amount of r-process material ejected per supernova, 
~ 10~ 5 M , as expected from the supernova rate and total r-process yields of the Galaxy. 
Finally, the conditions in the wind are determined by the properties of the neutron star and 
the characteristics of its neutrino emission. Therefore, the r-process nucleosynthesis might 
be approximately constant from event to event, reflecting the near constancy of the neutron 
star mass and cooling history, if possible complications due to fallback could be ignored. 

However, problems have emerged in the theoretical model. The successful r-process 
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calculations, with the possible exception of Woosley et al. (1994), all utilize parametric 
modifications to the key parameters, especially an artificial increase in the entropy. The 
high entropies of Wilson's supernova model reported in Woosley et al. (1994) have not been 
replicated elsewhere. This discrepancy has been emphasized by the analytic calculations of 
Qian & Woosley (1996, hereafter Paper I). In general, the deficiency in entropy is only a 
factor of two, but the gap is proving difficult to bridge. Paper I also suggested that other 
relevant parameters of the problem, especially a short dynamic time scale or large neutron 
excess, might compensate for the low entropy, or that there may emerge other sources of 
entropy hitherto neglected. But, for the time being, a potentially beautiful solution to a 
classic problem — the origin of the r-process — falters by a factor of two. 

We will not resolve this quandary in the present paper. What we shall do however, 
is: a) show that the standard models of the neutrino-driven wind derived in Paper I give 
interesting nucleosynthesis above the iron group, but with the electron fractions obtained 
in Woosley et al. (1994), do not give the classical r-process; and b) determine, both 
analytically and numerically, the conditions that are required in this sort of model to 
produce the r-process, in particular the platinum peak. 

2. Nucleosynthesis In Neutrino-Driven Winds 

An r-process might occur for various combinations of entropy, electron fraction, and 
dynamic time scale. Different authors, using numerical supernova models, have arrived at 
qualitatively different values for these key parameters in the neutrino-driven wind. In order 
to obtain a better understanding of the physical conditions, both analytic and numerical 
studies of the wind were carried out in Paper I. The primary goal of that paper was to 
examine the dependencies of the physical parameters in the wind on the neutron star mass 
(M) and radius (R), and the emergent neutrino luminosity (L u ) and energy spectra. The 
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analytic study was in the same spirit as that of Duncan, Shapiro, & Wasserman (1986), 
but was more extensive and directed towards nucleosynthetic issues. The analytic results of 
Paper I were given by equations (48a), (48b), and (49) for the entropy (S), equations (58a) 
and (58b) for the mass outflow rate (M), equation (61) for the dynamic time scale (Td yn ), 
and equation (77) for the electron fraction (Y e ). 

Physically, the first three parameters (S, M, and Td yn ) are determined by the sum of 
heating produced by all neutrinos, whereas the evolution of Y e mostly reflects the difference 
between the v e and v e fluxes through the inter-conversion of free nucleons by v e and v e 
captures. As discussed in Paper I, the first three parameters are not sensitive to the exact 
values of Y e , and their determination can be essentially decoupled from the evolution of 
Y e in the wind. The analytic results for these three parameters were tested by a series of 
numerical calculations using the one-dimensional implicit hydrodynamic code KEPLER. 
A total of nine models were studied. The neutron star mass and radius, and neutrino 
luminosity at the inner boundary of these models were varied to observe the corresponding 
effects on these three parameters. A comparison between analytic and numerical results 
was summarized in Table 1 of Paper I. 

From these numerical models, we have extracted the velocity (v), density (p), and 
temperature (T) of a mass element as functions of its position (r) in a steady-state wind. 
Using dt = dr/v(r), we can obtain the evolution of r, p, and T with time t for the mass 
element. Starting at T 9 « 10 (T 9 is the temperature in units of 10 9 K), when nuclear 
statistical equilibrium (NSE) is assured, we then follow the time evolution of the nuclear 
composition in this mass element with a reaction network ( |Woosley fc Hoffman 1992|) , until 
the composition freezes out at Tg ~ 1. To good accuracy, the mass element is initially 
composed of free nucleons, in proportions specified by the initial electron fraction Y e ti . The 
initial electron fraction is essentially determined by the equilibrium between v e and P e 
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captures on free neutrons and protons, respectively ( |Qian et al. 1993| ; see also Paper I). If 



we denote the rates for v e and v e captures on free nucleons as X Uen and A PeP , respectively, 
the initial electron fraction is given by 

In turn, the rate X UeU (Ap eP ) is determined by the v e (u e ) luminosity and energy spectrum. 

In realistic supernova models, the neutrino luminosity and energy spectra evolve with 
time. The individual wind models of Paper I represent the steady-state configurations 
reached at different neutrino luminosities over time scales much shorter than the evolution 
time scales of the neutrino luminosity and energy spectra. Because the time evolution of 
the neutrino energy spectra is much less pronounced than that of the neutrino luminosity, 
two generic sets of neutrino mean energies were assumed for these models. From the 
analytic results of Paper I, we can see that the mass loss rate, the dynamic time scale, and 
especially the entropy would not change significantly if a more precise prescription of the 
neutrino energy spectra were used. However, because the electron fraction is sensitive to the 
difference between the v e and v e fluxes, and the nature of heavy element nucleosynthesis is 
extremely sensitive to Y e , in this paper we calculate the v e and v e reaction rates according to 
the time evolution of neutrino luminosity and energy spectra in Wilson's 20 M supernova 
model used in Woosley et al. (1994). Specifically, we take from Wilson's supernova model 
the neutrino energy spectra corresponding to the same neutrino luminosity as in the wind 
model. Starting with the initial value in equation (1), we then follow the evolution of Y e in 
the reaction network, taking into account u e and v e captures on free nucleons and heavy 
nuclei ( |McLaughlin Fuller 1991) ), electron and positron captures on free nucleons and 



heavy nuclei, and nuclear /3-decays. 

The results of the nucleosynthesis calculations are presented in Tables 1-3. The 
entropy, the initial electron fraction, and the dynamic time scale are given for each wind 
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model. The dynamic time scale roughly corresponds to the time over which the temperature 
changes by one e-fold (Paper I). In order to check the influence of the neutrino flux on 
the nucleosynthesis, we have carried out five different calculations. All runs included 
electron and positron captures on free nucleons and nuclei, as well as nuclear /3-decays. 
The individual runs differ in the inclusion of various neutrino reactions. Respectively, they 
cover the cases including (1) no neutrino reactions (column 2), (2) v e and v e captures on 
free nucleons only (column 3), (3) v e and 9 e captures on free nucleons, and neutral-current 
neutrino spallation on a-particles (column 4), (4) v e and v e captures on free nucleons 
and nuclei (column 5), and (5) u e and v e captures on free nucleons and nuclei, as well as 
neutral-current neutrino spallation on a-particles (column 6). For all runs, we give the 
electron fraction {Y e j), the average mass number of nuclei excluding free nucleons and 
a-particles (A), and the neutron and a-particle mass fractions (X n j and X a j) at the 
freeze-out of the charged-particle reactions (Tg ~ 2.5). 

In Figures 1-9, detailed nucleosynthesis results from the runs that included u e 
and v e captures on free nucleons (column 3 in Tables 1-3) are given in terms of the 
production factor, defined as the final mass fraction (after all weak decays) of a given stable 
nucleus divided by its solar abundance flAnders fc Grevesse 1989j ). These results are also 



representative of the nucleosynthesis obtained in the other runs. In these figures, the most 
abundant isotope in the solar abundance distribution for a given element is plotted as an 
asterisk. Isotopes of a given element are connected by solid lines. A diamond around a 
data point indicates that the isotope is produced chiefly as a (neutron-rich) radioactive 
progenitor. The dotted horizontal lines represent an approximate "normalization band," 
bounded from above by the largest production factor in the calculation and from below by a 
production factor four times smaller. Nuclei that fall within this band will be the dominant 
species produced. In Figures 1-9, no re- normalization has been attempted. 
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Models 10A-F produced interesting nucleosynthesis representative of the a-process 
( |Woosley fc Hoffman 1992| ; [Witti, Janka, fc Takahashi 1994|) . The most abundant nuclei 
produced have mass numbers 90 < A < 110. Model 10A (Y e « 0.47, S ~ 70 per baryon) 
shows the production of the N = 50 closed-neutron-shell nuclei which were grossly 
overproduced in previous studies. Models 10B and C, with progressively lower neutrino 
luminosity and lower values of Y e , made heavier nuclei. Production of Sn, Sb, and Te was 
not accurately calculated, as the radioactive progenitors for these species were isotopes of 
Ru, the last element in our reaction network. Models 30A-C produced nuclei near the iron 
group. Model 30A exhibits interesting nucleosynthesis for Y e > 0.5, while Model 30C shows 
the production of 64 Zn (made as itself), the dominant isotope of this element. This nucleus 
was not accounted for in the surveys of Galactic chemical evolution and nucleosynthesis in 
massive stars QTimmes, Woosley, fc Weaver 199"5"| |Woosley fc Weaver 1995| ), and appears to 



be made predominantly under conditions similar to those obtained in the neutrino-driven 
wind ( Hoffman et al. 19961 ). 



From Tables 1-3, it is also clear that inclusion of the neutrino reactions made a 
difference. The v e and v e captures on free nucleons have the largest effect. The electron 
fraction increases appreciably due to these capture reactions when free nucleons are being 
assembled into a-particles. This so-called "a-effect" (|Fuller fc Meyer 1995|; |McLaughlin 



[Fuller, fc Wilson 1996| ) is evident when we compare the case including no neutrino reactions 
(column 2) with the cases including various neutrino reactions (columns 3-6). The inclusion 
of neutral-current neutrino spallation on a-particles and v e and v e captures on heavy 
nuclei did not have a major effect on the nucleosynthesis, at least before the freeze-out of 
the charged-particle reactions. For the relatively low entropies studied here, the neutrino 
spallation on a-particles did not have any appreciable influence on the final a-particle mass 
fraction. This is to be contrasted with the dramatic effect of these spallation reactions on 
the r-process in Wilson's high-entropy supernova model QMeyer 1995| ). The v e captures on 
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heavy nuclei may have important consequences for the ensuing neutron-capture phase of 
the r- process ( [Nadyozhin fc Panov 1993|) . Regardless of these issues, it is clear that none of 



these wind models produced an r-process. The neutrons had been consumed by the time 
the charged-particle reactions froze out (Tg « 2.5), and the nuclear flow did not even reach 
the A ~ 130 r-process peak. 

Another important issue concerns the overall ejection of the synthesized material in 
the wind into the interstellar medium. As the following simple argument will show, the 
magnitudes of the production factors may preclude the ejection of such material in all nine 
wind models. From the time evolution of neutrino luminosity in Wilson's supernova model, 
we find that, for example, the neutrino luminosity decreases from twice to half the value 
in Model 30B over a time of r u ~ 1 s. With a mass loss rate of M ~ 1.1 x 10~ 2 M & s _1 
and the largest production factor of (X w /X Q ) max rs 2.7 x 10 4 in Model 30B (Figure 8), the 
corresponding normalized production factor is ~ (Mr^/20 M Q )(X w /X Q ) mauX ~ 15, if the 
total amount of ejecta from the supernova is 20 M & . Models 30A and C give normalized 
production factors of up to ~ 30. For Models 10A-F, the normalized production factors for 
the nuclei produced in the largest amount are typically of order 100. However, Woosley & 
Weaver (1993) find that the normalized production factor should not be much above 10 in 
order for supernovae to produce the observed solar abundance of oxygen. Therefore, these 
wind models cannot represent in detail what commonly occurs in supernovae. 

It is quite possible that appropriate modifications of these standard wind models can 
lead to the physical conditions for acceptable nucleosynthesis. Paper I studied the effects of 
an additional energy source on the entropy and dynamic time scale in the wind. Model 10F 
of Paper I was recalculated with an additional energy input of 5 x 10 47 erg s _1 distributed 
uniformly in volume between 15 and 25 km. At these radii, the mass loss rate has been 
more or less determined. The additional energy input represents a moderate perturbation 
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to the total amount of heating provided by neutrinos (1.2 x 10 48 erg s _1 ). As a result, the 
entropy increased from 140 to 192 per baryon, and the dynamic time scale decreased from 
0.11 to 0.022 s, while the mass loss rate slightly increased from 2.8 x 10~ 6 to 3.7 x 10~ 6 M & 
s" 1 . 

The effect on the nucleosynthesis was dramatic as shown in Table 4. The neutron-to- 
seed ratio (cf. eq. [3]), less than 10 in all of the unmodified models, rose to ~ 166. With an 
average mass number of A ~ 90 for the seed nuclei, uranium could be produced if all the 
neutrons were to be subsequently captured. The A ~ 195 r-process peak probably could 
have been produced with less additional energy input, and hence a lower entropy and a 
longer dynamic time scale, than assumed in Paper I. We conclude that lower values of Y e 
than those calculated by Wilson in Woosley et al. (1994), or additional energy input like 
that considered in Paper I, are necessary to produce an r-process in a spherically symmetric 
wind model. 



3. The Requisite Conditions For The r-Process 

From our studies of nucleosynthesis in the neutrino-driven wind in the previous section, 
we learn that the standard wind models of Paper I fail to make an r-process. On the other 
hand, we also see that reasonable modifications of these models can significantly change 
the physical conditions in the wind, and therefore give rise to a possible r-process. In order 
to better understand the deficiencies of the standard wind models, and furthermore, to 
motivate and direct physically plausible modifications of these models, we now survey the 
important physical parameters required for a strong r-process. (See also |Takahashi, Witt 



fe Janka 1991 ; ITakahashi 1996j ; Freiburghaus et al. 19961 ; and |Meyer fc Brown 19961 ). 



We consider the following generic model for the r-process. Neutron-rich material 
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initially composed of free nucleons at high temperatures (T 9 9) adiabatically expands and 
cools. After nearly all the protons are assembled into a-particles at T 9 ss 5, an a-process 
occurs to burn the a-particles into heavy nuclei. The a-process stops when charged-particle 
reactions freeze out at T 9 ps 2.5. The heavy nuclei produced at the end of the a-process 
then become the seed nuclei for the subsequent rapid neutron capture process, or the 
r-process. We do not intend to account for the full detail of the r-process, such as the 
final abundance distribution, which becomes meaningful only in the context of a consistent 
astrophysical model. What we are most interested in is the physical conditions favorable 
for the production of the most abundant r-process nuclei, such as those in the platinum 
peak of the solar r-process abundance distribution. For this purpose, we can think of the 
r-process as the transformation of seed nuclei into r-process nuclei through simple addition 
of the available neutrons. In this sense, the possibility of producing r-process nuclei around 
a certain mass number depends only on the relative abundances of seed nuclei and neutrons 
at the end of the a-process. 

In general, the composition resulting from the a-process satisfies 

X nJ + X aJ + X a & 1, (2) 

where X n j and X a j are the final mass fractions of neutrons and a-particles, respectively, 
and X s is the total mass fraction of seed nuclei. If we represent the seed nuclei with a mean 
proton number Z and a mean mass number A, we can define a neutron-to-seed ratio at the 
end of the a-process as 

n ~ XnJ A. (3) 



s X s 

In terms of this ratio, our simplified condition for making r-process nuclei with mass 
number A becomes 

Tl — 

- + A«A (4) 

s 
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We are particularly interested in the production of the platinum peak, and will describe the 
numerical calculations for A ps 200 in the following. 

To follow the nucleosynthesis in the adiabatically cooling material with a nuclear 
reaction network, we need the initial composition and the time evolution of temperature 
and density. The initial composition at T 9 ps 9 can be simply specified by the initial 
electron fraction Y €ti , with the initial mass fractions of neutrons and protons given by 
X n>i wl - Y e>i and X Pji ps Y e>i , respectively. Because temperature and density are related 
through the constant entropy for the adiabatically cooling material, we only need to specify 
the temperature as a function of time. For simplicity, we introduce a dynamic time scale 
(Tdyn) over which the temperature changes by one e-fold, i.e., 

T 9 (t)«T 9 (0)exp(-t/r dyn ). (5) 

With this time evolution of the temperature, the duration of the adiabatic expansion from 
T 9 w 9 to 2.5 is 

t cxp ss 1.28r dyn . (6) 

Hereafter, we will refer to t cxp as the expansion time. We further assume that the entropy 
is dominated by contributions from radiation and relativistic electron-positron pairs, and is 
given by 

5^3.33^, (7) 
Ps 

where p 5 is the density in units of 10 5 g cm -3 , and S is in units of Boltzmann constant per 
baryon. For convenience, we frequently refer to entropy without its unit throughout this 
paper. 

Now it is straightforward to determine the combinations of initial electron fraction, 
entropy, and expansion time, for which an a-process can lead to a sufficient neutron-to-seed 
ratio for production of r-process nuclei with A ps 200. We choose a range of expansion 
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times (0.005 < t cxp < 0.25 s). For each t exp , we survey a broad range of initial electron 
fractions {Y e>i = 0.20, 0.25, 0.30, 0.35, 0.40, 0.45, 0.46, 0.47, 0.48, 0.49, and 0.495). With 
a particular set of t exp and Y Si i, we seek through iteration the appropriate entropy which 
enables the a-process to produce a final composition satisfying equation @ for A « 200. 
In our calculations, we take into account the effects of electron and positron captures and 
nuclear /3-decays on the evolution of the electron fraction. 

The results are given in Table 5 and Figure 10. For each successful run, the 
composition at the end of the a-process are given in terms of the final neutron and 
a-particle mass fractions and the mean proton and mass numbers for seed nuclei. The 
mean proton and mass numbers were calculated from Z = Y(Z S , A S )Z S / Y(Z S , A s ) and 
A = Y(Z S , A S )A S / Y(Z S , A s ), respectively, with Y(Z S ,A S ) the number fraction of the 
seed nucleus with proton number Z s and mass number A s . The final electron fraction and 
the neutron-to-seed ratio are also given in Table 5. The combinations of initial electron 
fraction and entropy in the successful runs for three specific expansion times are shown as 
filled circles connected by solid lines in Figure 10. At a given t exp , values of Y e ^ and S to 
the left of the solid line will not give a sufficient neutron-to-seed ratio for production of the 
platinum peak. In addition, the following features of this figure are worth mentioning: 

(1) Depending on the expansion time, there exist many possible combinations of 
initial electron fraction and entropy that can produce nuclei with A ~ 195. The "high 
entropy" scenario (S ^ 350 and 0.495 ^ Y e j ^ 0.40) corresponds to longer expansion times 
(t cxp ^ 0.1 s). The results for t exp = 0.25 s are consistent with the conditions seen at 
late times in Wilson's supernova model flWoosley et al. 1994| ), and those employed in the 
successful r-process calculations of Woosley et al. (1994) and Takahashi, Witti, & Janka 
(1994), although the later required an artificial increase in the entropy (by a factor of five) 
to produce nuclei with A ~ 195. 
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(2) Alternatively, there is a "low entropy" scenario (S <; 200 and Y Cii <^ 0.4) that 
requires shorter expansion times (t cxp <; 0.025 s). For the most extreme case shown 
(^cx P = 0.005 s), the platinum peak can be made for any Y e>i ^ 0.495 if S ~ 150. Such 
an expansion would correspond to very high velocities which might be more appropriate 
to a jet than to a quasi-steady-state wind as considered in Paper I. However, these rapid 
expansions may not continue after the a-process in order to allow enough time for the 
neutron capture phase of the r-process. The slowing down of the expansion could be 
facilitated by a massive overlying mantle in the case of a Type II supernova. 

(3) For smaller values of Y Cii , the required entropy decreases regardless of the expansion 
time. This merely reflects the very neutron-rich nature of the initial composition. The 
results of Paper I suggest that values of Y e below 0.3 might be very difficult to achieve in 
the neutrino-driven wind. Material with Y e < 0.3 would need to be ejected without any 
significant interaction with neutrinos. 

(4) For a fixed expansion time, the required value of entropy actually decreases as Y e>i 
increases from ~ 0.48 to ~ 0.495. Lower entropy translates to higher density (cf. eq. [7]), 
and would normally produce more seed nuclei. However, in these cases of high Y" e i , the 
electron fraction has even more leverage on the seed production. As Y e ^ increases towards 
0.5, the neutron abundance decreases to vanishingly small values when the a-process begins 
at T 9 pa 5. This in turn diminishes the efficiency of burning a-particles through the main 
reaction path bridging the unstable mass gaps at A = 5 and 8, i.e., 4 He(o;n, 7) 9 Be(o;, n) 12 C. 
Table 5 shows that these high Y e>i values give extreme a-rich freeze-outs with final a-particle 
mass fractions approaching unity and comparable mass fractions for neutrons and heavy 
seed nuclei. 

In deriving the above results, we have made two major assumptions: (1) the entropy is 
proportional to T 3 / p with the proportionality constant calculated for a mixture of radiation 
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and relativistic electron-positron pairs, and (2) weak interactions other than electron and 
positron captures and nuclear /3-decays can be neglected in the nuclear reaction network. 
We will discuss how our results are affected if we drop either assumption in turn. 

Due to its logarithmic dependence on the temperature and density, the entropy 
of non-relativistic particles roughly stays constant over the temperature range in our 
calculations (9 ^ Tg ^ 2.5). For a total entropy of S ^ 20, the change in density 
essentially maintains a constant entropy of relativistic particles as the material adiabatically 
cools. At high temperatures (9 ^ Tg ^ 5), the relativistic particles include photons and 
electron-positron pairs, and the corresponding entropy is given by equation (7). However, 
as the temperature cools below Tg ~ 5, electron-positron pairs begin to annihilate. The 
situation is much like the Big Bang. Eventually, at Tg ~ 1, the only relativistic particles 
in the material are photons, with the corresponding entropy given by S ~ 1.21Tg/p 5 . In 
general, we can write the entropy in relativistic particles as S ~ C(Tg)Tg / p 5 . Because 
C(Tg) decreases noticeably from Tg « 5 to 2.5 when the a-process is taking place in 
our calculations, we have overestimated the density by using C{Tg) « 3.33 throughout 
the adiabatic expansion of the material. Consequently, we have overestimated the seed 
production and underestimated the neutron-to-seed ratio at the end of the a-process. This 
is especially true for the cases of high entropies where more time is available to produce 
seed nuclei. 

With an accurate expression for the entropy, the condition in equation (^) is satisfied at 
slightly lower entropies than those given in Table 5 for the same initial electron fraction and 
expansion time. In fact, we have repeated our calculations for Y e ^ = 0.30, 0.40, 0.45, 0.47, 
and 0.49, using an exact adiabatic equation of state to compute the density corresponding 
to a specific temperature. This equation of state takes into account the contributions to 
the entropy from radiation, electron-positron pairs, and ions. The results are shown as 
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filled squares in Figure 10 for three of the expansion times explored in the initial numerical 
survey. As expected, our previous results using equation (7) overestimated the entropy 
required to produce the platinum peak. When we use the exact adiabatic equation of state, 
the required entropy is lower by ~ 10% for the longest expansion time t exp = 0.25 s, whereas 
for the shortest expansion time t cxp = 0.005 s, the results are essentially unchanged. 

If the r-process occurs in an environment with intense neutrino flux, perhaps we 
should also include various neutrino interactions in the nuclear reaction network. In fact, 
we could have specified a less "generic" model by considering an r-process site similar 
to the neutrino-driven wind. In this case, the adiabatically expanding material is also 
moving away from the neutrino source. We can define a constant dynamic time scale as 
r d y n ~ f/v, with r the distance from the neutrino source and v the expansion velocity. 
The time evolution of temperature in equation (5) follows on assuming T oc r _1 . We can 
then introduce an additional parameter, e.g., the initial neutrino flux $ ui at T 9 « 9, in our 
calculations. As the material adiabatically expands, the neutrino flux it receives decreases 
as r~ 2 oc exp(— 2t/rd yn ). In principle, using the above prescription, we can repeat our 
calculations for a range of Q vi with various neutrino interactions included in the nuclear 
reaction network. 

From our discussions in §2, we have seen that for relatively low entropies of S <^ 200, 
the major role of the neutrino flux in determining the neutron-to-seed ratio is to increase 
the electron fraction by u e and V e captures on free nucleons through the ct-effect. In 
addition, Meyer (1995) has shown that for high entropies of S ~ 400, neutral-current 
neutrino spallations on a-particles during the a-process can increase the production of seed 
nuclei. In both cases, inclusion of neutrino interactions tends to reduce the neutron-to-seed 
ratio at the end of the a-process, although the effects of these neutrino interactions are 
less important for shorter expansion times. Therefore, we can interpret the entropies in 
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Table 5 and Figure 10 as the minimum values required to produce the platinum peak for 
given sets of initial electron fraction and expansion time. With this interpretation, we can 
avoid repeating our calculations and complicating our results with an additional parameter 
Of course, in a consistent astrophysical model for the r-process where intense neutrino 
flux exists, the exact conditions for production of, e.g., the platinum peak, have to be 
determined with full consideration of various neutrino interactions. 



4. Analytic Treatment Of The a-Process 

In order to provide some physical insight into what determines the neutron-to-seed 
ratio, and extend our results in the previous section to production of r-process nuclei in 
general, we present an analytic treatment of the a-process in this section based on our 
generic r-process model. If we ignore possible neutrino interactions, the electron fraction at 
the beginning of the a-process (T 9 5) is about the same as the initial electron fraction 
F e ,j at T 9 pa 9. At T 9 ~ 5, the material is essentially composed of free neutrons and 
a-particles for Y Cii < 0.5, with almost all the protons already assembled into the a-particles. 
The mass fractions of a-particles and neutrons at the beginning of the a-process are then 
approximately given by 

X a>0 « 2Y eji , (8a) 
X nfi » l-2Y e>i , (8b) 

respectively. 

The composition at the end of the a-process (T 9 ps 2.5) satisfies 

l -X aJ + jX s » Y eJ , (9) 

where Y e j is the final electron fraction at T 9 ps 2.5. Using equation (9) together with 
equations (2)-(4), we find that the final mass fractions of a-particles and neutrons have to 
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be 

**f « ^f'^t (10a) 
a,/ 1/2 -Z/A' v ; 

*n,/ « ^^ (1-21;,/), (10b) 

respectively, in order to produce r-process nuclei with mass number A. Comparing 
equations (8a) and (8b) with equations (10a) and (10b), we see that the fractional change 
in the neutron mass fraction during the a-process is less than that in the a-particle mass 
fraction for Y e j pa Y e j < Z/A. Obviously, because neutrons carry no charge and the 
mean charge per nucleon for a-particles exceeds Z/A for the heavy seed nuclei, the final 
composition favors the presence of neutrons for Y e j < Z/A. Accordingly, we will present 
the analytic treatment of the a-process for two different cases: Y £ti < Z / A and Y £ti > Z/A. 
In both cases, we will assume Y e j pa Y €ti . 

As the temperature declines from T 9 pa 5 to 2.5, a-particles and neutrons are partially 
assembled into heavy seed nuclei. We can describe the a-process in terms of the time 
evolution of the a-particle and neutron abundances. During the a-process, the burning 
of a-particles mainly proceeds via the reaction sequence 4 He(an, 7) 9 Be(a, n) 12 C. The 
production of seed nuclei occurs through the efficient a-capture reactions starting with 
9 Be(a,n) 12 C. Consequently, the rates of change in the a-particle and neutron number 
fractions can be approximately written as 

dY 

pa -FY a Y 9 pN A (av) an , (11a) 

dY 

Pa -GY a Y 9 pN A (av) an , (lib) 

respectively, where Y 9 is the number fraction of 9 Be, and NA{<rv) an is the reaction rate 
for 9 Be(a,n) 12 C in units of cm 3 s _1 g -1 . In equations (11a) and (lib), F and G are the 
numbers of a-particles and neutrons that make up a typical heavy seed nucleus. For a seed 
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distribution with mean proton number Z and mean mass number A, we have F ps Z/2 and 
A-2Z. 

Due to its low Q-value of only 1.573 MeV, the reaction 4 He(o;n, 7) 9 Be is tightly- 
balanced by its reverse reaction essentially over the entire temperature range 5 ^ T 9 ^ 2.5. 
According to statistical equilibrium, the number fraction of 9 Be during the a-process is 
given by 

Y 9 « 8.66 x 10- n F a 2 y;P5 T 9" 3 exp(18.26/T 9 ). (12) 

Because the density p is related to the temperature T 9 through the constant entropy S, 
and NA{(Jv) an depends on T 9 only, equations (11a) and (lib) can be expressed in a more 
convenient form if we regard Y a and Y n as functions of temperature. Using equation (5), we 
have 

dY 

" ~ FY^Y n g(T 9 )r dyQ , (13a) 



dT 9 
dY n 



GY^Y n g(T 9 )r dyn , (13b) 



dT 9 

where g(T 9 ) has the unit of s _1 , and is given by 

g(T 9 ) « 8.66 x 10~ 6 p 3 5 T 9 4 exp(18.26/T 9 )N A (av) an . (14) 

Now we can solve equations (13a) and (13b) for the two different cases mentioned 
previously. For Y e>i < Z / A, the final composition favors the presence of neutrons. So we can 
approximately take Y n m Y n $ = X n>0 during the a-process. In this case, it is straightforward 
to solve equation (13a) and obtain 

Y-J - Yj » 2FY nfi r dyn £ g(T 9 )dT 9 . (15) 

Likewise, for Y e>i > Z / A, we can approximately take Y a ss Y a $ = X a $/A during the 
a-process, and solve equation (13b) to obtain 



Y n j ~ Y nfi exp 



-GYl r dyn f g{T 9 )dT 9 

J2.5 



(16) 
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Equations (15) and (16) implicitly constrain the combinations of Y £ti , S, and 7d yn , for which 
the a-process can give a sufficient neutron-to-seed ratio for production of r-process nuclei 
with mass number A. 

To proceed further, we approximate the constant entropy during the adiabatic 
expansion as S ~ C(T 9 )Tg / p 5 , with C{Tg) decreasing from 3.33 at T 9 ^ 5 to 1.21 at T 9 <; 1. 
So equation (14) can be rewritten as 

g(T g ) w 8.66 x 10~ 6 S- 3 C(T 9 )*I$ exp(18.2Q/Tg)N A (av) an . (17) 

We plot g(T 9 )S 3 as a function of T 9 in Figure 11, assuming C(T 9 ) pa 3.33 and using the 
fitting formula for NA(<Jv) an given by Wrean, Brune, & Kavanagh (1994). As we can 
see from this figure, g(Tg) decreases monotonically with temperature over 2.5 <^ T 9 <; 5. 
At T 9 w 4, g(T 9 ) has already fallen below half its value at T 9 pa 5. Clearly, the main 
contribution to the integral J 2 5 5 g(T 9 )dTg comes from 4 ^ T 9 <^ 5. This remains true even 
when the exact form of C(Tg) is used. For our analytic estimates, we can approximately 
evaluate the integral with C(T 9 ) pa 3.33, and obtain 

/ 5 g{Tg)dTg pa 6.4 x 10 8 S~ 3 s" 1 . (18) 

J2.5 

Using equation ( p!8|) and assuming F e j pa Y" e i; we can rearrange equations (15) and (16) as 

4 x 10 7 Z(1 - 2Y eA 



1/3 



5 



5 



'X/2-Z/A' 


2 


1 


[Y etl -Z/A_ 




. 25^^ _ 



Tdyn 



for y eji < -=, 



VI. 



8 x 10 7 (A-2Z) /r dyn 



1/3 



In 



1-2Z/A 
1-A/A 



for y e>i > 



(19a) 
(19b) 



where we have also used equations (8a), (8b), (10a) and (10b). 

To compare our analytic results in equations (19a) and (19b) with the numerical results 



for A pa 200 in §3, we take Z pa 34 and A pa 90 from the numerical survey, and obtain 

1/3 



s 



10 d 



1 - 


2V^j 






0.33 


2 


l 


2 


L(y e , 8 -o.i7)_ 









cxp 



, for Y e4 < 0.38, 



(20a) 
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( 



t, 



) 



1/3 



s 



2 x io 3 y M 



exp 



S 



for Y^ t > 0.38, 



(20b) 



where we have replaced Td yn with t cxp « 1.28rd yn . These analytic results are shown as open 
circles connected by dotted lines in Figure 10. The general agreement between our analytic 
results and the numerical survey using equation (7) (filled circles connected by solid lines) 
is quite good. Because Y a and Y n can decrease during the a-process and effective burning 
of a-particles may start at T 9 < 5, we always tend to overestimate the entropy by holding 
either Y a or Y n constant and doing the integral in equation ( Jl8| ) over the entire temperature 
range 2.5 <^ Tg ^ 5 in our analytic treatment. The largest discrepancies occur at Y e ^ 0.47 
where we overestimate the entropy by ~ (15-50)%. Using the specific values of Z and A 
found in the numerical survey would slightly improve the agreement. Further improvement 
could be obtained by solving equations (11a) and (lib) together instead of approximately 
solving each for a specific case. However, these improvements would add little to our 
understanding of the physics determining the neutron-to-seed ratio. We also notice that the 
same level of agreement with the numerical survey using an exact equation of state (filled 
squares) can be achieved if we use C(Tg) ~ 3 instead of 3.33 to account for the annihilation 
of electron-positron pairs into photons at T 9 < 5. 

From our analytic treatment of the a-process, we can clearly see the individual 
roles of the initial electron fraction, entropy, and dynamic time scale in determining the 
neutron-to-seed ratio. In addition to specifying the overall availability of neutrons (cf. 
eq. [8b]), the initial electron fraction serves to direct the path of nuclear flow during 
the a-process. As our analytic treatment indicates, the comparison of Y e i with the ratio 
Z/A for the typical seed distribution reflects whether a-particles or neutrons are mainly 
consumed during the a-process. The influence of entropy is manifested through the density 
dependences of the equilibrium abundance of 9 Be (cf. eq. [12]) and the rate for burning 
a-particles (cf. eq. [11a]). Physically, a high entropy means many photons per baryon 
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in radiation-dominated conditions. A significant fraction of these photons can be on the 
high-energy tail of the Bose-Einstein distribution, and therefore can maintain a low 9 Be 
abundance through the photo-disintegration reactions. In turn, this limits the overall 
efficiency of burning a-particles, and hence the production of seed nuclei. The dynamic time 
scale, or the expansion time, specifies the duration of the a-process (cf. eq. [6]). Obviously, 
the expansion time also acts to limit the production of seed nuclei. In general, a lower 
initial electron fraction, a higher entropy, and a shorter expansion time all give a larger 
neutron-to-seed ratio. It is most interesting to notice that for a given Y £ti , the composition 
resulting from the a-process essentially only depends on the combination S 3 /t cxp (cf. eqs. 
[20a] and [20b]). Consequently, the same neutron-to-seed ratio can be achieved for the same 
initial electron fraction with an expansion time 8 times shorter if the entropy is reduced by 
a factor of 2. 



5. Conclusions 

For reasonable assumptions regarding neutrino luminosity, neutron star mass and 
radius, and the time history for Y e , the nucleosynthesis resulting from the analytic model 
developed in Paper I for the neutrino-driven wind does not resemble the solar r-process, 
although a number of interesting species in the mass range 90 < A < 120 are produced. 
This failure may be a consequence of important (but unknown) physics neglected in Paper 
I, or our results may reflect the true nucleosynthesis from typical core-collapse-driven 
supernovae (however, see the discussion concerning the ejection of the wind material in §3). 
Extra (but, so far, artificial) energy input to the wind beyond the injection radius (where 
M is determined) does give a successful r-process. Possible sources of this energy were 
discussed in Paper I. 

A numerical survey has delineated the necessary combinations of key parameters 
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- Y e , entropy, and expansion time — needed to produce the third r-process peak (i.e., 
the platinum peak). High entropy is not a unique requirement. A shorter expansion 
time also serves to limit the number of heavy seed nuclei produced and thus increase the 
neutron-to-seed ratio. A lower Y e also leads to a larger neutron-to-seed ratio. The sensitivity 
of the neutron-to-seed ratio to Y e diminishes as one proceeds to shorter expansion times, as 
does the sensitivity to the entropy. Specific values of Y e , entropy, and expansion time to 
produce the third r-process peak are given in Table 5 and Figure 10. 

Approximate analytic formulae (eqs. [20a] and [20b]) were derived that give the 
requisite entropy needed to produce the heavy r-process nuclei as functions of Y e and 
expansion time. These equations can be used to gauge whether other unstudied supernova 
models or other astrophysical environments are appropriate sites for making r-process 
nuclei. 

Given that the standard wind models, without artificial modification, produce a set 
of abundances distinct from the r-process, one must be concerned about the observational 
consequences. One possibility already mentioned is that important physics has still been 
omitted from the simple wind model — e.g., neutrino flavor mixing, added energy input 
from shocks, rotation, or magnetic fields, convection, etc. - and that the conditions 
required for the r-process may still ultimately be achieved in common core-collapse-driven 
supernovae. Perhaps material having a much lower Y e than calculated by Wilson in Woosley 
et al. (1994) is ejected ( [Burrows, Hayes, fc Fryxell 1995| ). Another possibility which must 
be seriously considered however, is that the r-process has more than one important site and 
neutron exposure. 

Sneden et al. (1996) and Cowan et al. (1996) have observed elements attributed to 
the third r-process peak in the metal poor halo giant stars HD 126238 ([Fe/H]= —1.7) 
and CS 22892-052 ([Fe/H]= -3.1). In HD 126238, the scaled solar abundances of both 
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Os and Pt have been clearly observed. Both elements are made almost exclusively by the 
r-process and fit the solar r-process abundance pattern. Similar results hold for the more 
metal poor star CS 22892-052, although Pt was not observed and the detection of Os is less 
certain. Coupled with previous data, the fit to the solar r-process pattern for all elements 
between Ba and Os is striking, suggesting that, for this range of nuclei, the solar r-process 
abundance distribution appears to be made in its entirety in the progenitor(s) of these 
metal poor stars. This result argues for a primary production scenario. Due to the star's 
very low metallicity, especially so for CS 22892-052, the observed r-process abundance 
pattern probably arose from only a few supernovae. 

Interestingly, the abundance pattern for CS 22892-052 shows that elements in the 
first neutron capture peak (Sr, Y, and Zr) are below their scaled solar r-process fractions 
relative to Pt, Os, and Th, yet well above the the iron group. Thus locally the r-process 
abundances seem solar, but globally they are not. Sneden et al. (1996) suggest that the 
bulk of the solar abundance of Sr, Y, and Zr is due to the s-process, but these nuclei 
are also easily produced in the neutrino-driven wind (Figure 1 here; [Woosley fc Hoffman 



19921 ; |Witti, Janka, fc Takahashi 1994| ; |Woosley et al. 1994j) . This is consistent with the 



existence of two sources for the r-process, one responsible for production of the r-isotopes 
for Z < 56, including those in the N = 50 peak, and the other for the heavier elements, 
possibly operating in another site or a higher entropy version of the same site. 

A very different type of r-process would arise also for much shorter expansion times, as 
might occur in accretion induced collapse where the wind is not slowed down by a massive 
overlying mantle ( |Woosley fc Baron 1992| ; |Woosley fc Hoffman 1992] ) . A high entropy 
r-process with short expansion time could occur through ejection by relativistic jets in 
coalescing neutron stars ( [Ruffert et al. 19"96| ). With short expansion times, the duration of 



the neutron capture phase of the r-process may require special consideration. If the material 
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undergoing nucleosynthesis cannot be slowed down during the neutron capture phase, the 
r-process may have to be accelerated by v e captures on heavy nuclei (Nadyozhin & Panov 
1993). Ultimately observational signatures of these very distinct physical processes may be 
needed to resolve the nature and site(s) of the r-process. 
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Fig. 1.— Nucleosynthesis in KEPLER wind model 10A: M = 1.4 M , R = 10 km, 

Lutot = 1-8 x 10 52 erg s _1 . Nuclei with N = 50 closed neutron-shell dominate the 
nucleosynthesis. 

Fig. 2.— Nucleosynthesis in KEPLER wind model 10B: M = 1.4 M , R = 10 km, 
L utot = 6.0 x 10 51 erg s" 1 . 



Fig. 3.— Nucleosynthesis in KEPLER wind model 10C: M = 1.4 M , R = 10 km, 
L utot = 3.6 x 10 51 erg s" 1 . 



Fig. 4.— Nucleosynthesis in KEPLER wind model 10D: M = 2.0 M , R = 10 km, 
L utot = 1.8 x 10 52 erg s" 1 . 



Fig. 5.— Nucleosynthesis in KEPLER wind model 10E: M = 2.0 M , R = 10 km, 
L utot = 6.0 x 10 51 erg s" 1 . 



Fig. 6.— Nucleosynthesis in KEPLER wind model 10F: M = 2.0 M , R = 10 km, 
L utot = 3.6 x 10 51 erg s" 1 . 
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Fig. 7.— Nucleosynthesis in KEPLER wind model 30A: M = 1.4 M , R = 30 km, 

Lutot = 1-8 x 10 53 erg s~\ These conditions lead to the production of species in the iron 
group. 

Fig. 8.— Nucleosynthesis in KEPLER wind model 30B: M = 1.4 M , R = 30 km, 
L utot = 6.0 x 10 52 erg s" 1 . 

Fig. 9.— Nucleosynthesis in KEPLER wind model 30C: M = 1.4 M , R = 30 km, 

Lutot = 3.6 x 10 52 erg s -1 . The production of 64 Zn (made as itself), a unique signature 
of the neutrino-driven wind, suggests that such winds have likely occurred in nature. 

Fig. 10. — The combinations of Y e , entropy, and expansion time required for the production 
of the A ~ 195 r-process peak nuclei. Circles connected by lines are for various fixed 
expansion times. Shown are the values derived in the numerical study using equation (7) 
(filled circles) and those from the analytic approximation (eqs. [20a] and [20b], open circles). 
The filled squares represent results from the numerical survey that used an exact adiabatic 
equation of state. 

Fig. 11. — The temperature dependent function g(Tg) (eq. [17]). The main contribution to 
the integral of g(T 9 ) comes from T 9 > 4. 
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Table 1. Nucleosynthesis in Wind Models 10A-C 





no-u 


z/(np a ) 




z/(np, A Z c ) 


v(np,a, A Z) 








Model 10A: S ~ 


74, Y e:l = 0.465, r dyn 


= 0.024 s 






.465 


.470 


.470 


.470 


.470 


A 


85.1 


81.0 


80.7 


80.3 


80.7 


X a 


.456 


.501 


.506 


.508 


.504 


X n 


4.4(-12) d 


3.0(-12) 


3.4(-12) 


3.4(-12) 


3.3(-12) 






Model 10B: S ~ 


87, Y e , t = 0.372, r dyn 


= 0.066 s 






.373 


.393 


.395 


.396 


.394 


A 


98.5 


96.8 


96.6 


96.6 


96.8 


X a 


.167 


.200 


.204 


.205 


.201 


x n 


4.73(-2) 


1.42(-2) 


1.17(-2) 


1.07(-2) 


1.28(-2) 






Model 10C: S - 


- 94, Y e ,i = 0.354, r dyn 


= 0.11 s 






.355 


.386 


.390 


.391 


.387 


A 


99.8 


97.7 


97.3 


97.3 


97.7 


X a 


.135 


.173 


.178 


.180 


.174 




7.75(-2) 


2.10(-2) 


1.57(-2) 


1.42(-2) 


1.91(-2) 



a z/ e and & e captures on nucleons (np) 
b neutral-current neutrino spallations on a-particles (a) 
c u e captures on heavy nuclei ( A Z) 
d x{y) e i x 10 s 



Table 2. Nucleosynthesis in Wind Models 10D-F 



no-v z/(np a ) z/(np,o; b ) u(np, A Z c ) v(np,a, A Z) 







Model 10D: S ~ 


109, Y e ,i = 0.465, r dyn 


= 0.024 s 






.465 


.470 


.471 


.471 


.471 


A 


90.0 


89.0 


89.0 


88.9 


88.9 


X a 


.642 


.679 


.682 


.683 


.683 


X n 


1.78(-5) d 


1.92(-6) 


1.74(-6) 


1.22(-6) 


1.22(-6) 






Model 10E: S ~ 


129, Y e ,i = 0.372, r dyn = 


= 0.066 s 






.372 


.396 


.398 


.399 


.399 


A 


98.0 


97.2 


97.1 


97.2 


97.2 


X a 


.302 


.345 


.348 


.348 


.348 


x n 


9.54(-2) 


5.10(-2) 


4.80(-2) 


4.66(-2) 


4.66(-2) 






Model 10F: S - 


140, Y^ t = 0.354, r dyn 


= 0.11 s 






.355 


.390 


.393 


.393 


.393 


A 


99.1 


98.1 


97.9 


98.0 


98.0 


X a 


.262 


.315 


.320 


.321 


.321 


x n 


1.26(-2) 


5.69(-2) 


5.16(-2) 


4.99(-2) 


4.99(-2) 



a z/ e and & e captures on nucleons (np) 
b neutral-current neutrino spallations on a-particles (a) 
c u e captures on heavy nuclei ( A Z) 
d x{y) e i x 10 s 
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Table 3. Nucleosynthesis in Wind Models 30A-C 





no-v 


z/(np a ) 


z/(np,a b ) 


z/(np, A Z c ) 


u(np,a, A Z) 








Model 30A: S ~ 


24, Y^ t = 0.510, r dyn 


= 0.032 s 






.510 


.509 


.509 


.508 


.508 


A 


56.1 


56.1 


56.1 


56.1 


56.1 


X a 


.271 


.271 


.271 


.272 


.272 


X n 


3.(-18) d 


l.(-15) 


l.(-15) 


l.(-15) 


l.(-15) 






Model 30B: S ~ 


26, Y e , t = 0.488, r dyn 


= 0.075 s 






.488 


.489 


.489 


.490 


.489 


A 


58.1 


58.0 


57.9 


57.9 


58.0 


X a 


.220 


.223 


.224 


.225 


.223 


x n 


3.(-15) 


2.(-15) 


2.(-15) 


3.(-15) 


3.(-15) 






Model 30C: S - 


- 28, Y e ,i = 0.481, r dyn 


= 0.12 s 






.481 


.483 


.484 


.484 


.483 


A 


59.6 


59.1 


59.0 


58.9 


59.1 


X a 


.181 


.190 


.194 


.194 


.191 


x n 


7.(-15) 


7.(-15) 


7.(-15) 


9.(-15) 


l.(-14) 



a z/ e and & e captures on nucleons (np) 
b neutral-current neutrino spallations on a-particles (a) 
c u e captures on heavy nuclei ( A Z) 
d x{y) e i x 10 s 



Table 4. Nucleosynthesis in Modified Wind Model 10F with Extra Energy Input 



no-z/ z/(np a ) z/(np,o; b ) v(np, A Z c ) u(np,a, A Z) 





Modified Model 10F: 


S ~ 192, Y e , 


= 0.354, r dyn = 0.022 s 




YeJ 


.354 


.363 


.363 


.363 


.363 


A 


89.7 


91.0 


91.0 


91.0 


91.0 


X a 


.615 


.629 


.629 


.629 


.629 


x n 


.258 


.240 


.240 


.240 


.240 



a z/ e and V e captures on nucleons (np) 

b neutral-current neutrino spallations on a-particles (a) 

c u e captures on heavy nuclei ( A Z) 
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Table 5. 


Results from the Numerical Survey Using S r 


i 3.33T 9 3 /p 5 




^exp (s) 




S 




X nJ 


Z 


A 


n/s 


0.005 


0.200 


28 


0.074 


0.462 


35.8 


101.8 


101 




0.250 


58 


0.239 


0.396 


34.7 


97.0 


105 




0.300 


81 


0.391 


0.320 


33.9 


93.8 


104 




0.350 


103 


0.545 


0.242 


33.3 


91.3 


104 




0.400 


124 


0.696 


0.163 


32.9 


89.4 


103 




0.450 


145 


0.852 


0.083 


32.4 


87.2 


111 




0.460 


147 


0.882 


0.067 


32.2 


86.6 


112 




0.470 


144 


0.910 


0.050 


32.0 


85.8 


108 




0.480 


131 


0.942 


0.033 


32.2 


87.9 


116 




0.490 


117 


0.970 


0.016 


34.7 


97.1 


107 




0.495 


112 


0.983 


0.007 


41.1 


116.0 


85 


0.025 


0.201 


50 


0.081 


0.466 


36.1 


101.8 


105 




0.251 


95 


0.232 


0.393 


35.2 


97.9 


103 




0.301 


136 


0.393 


0.320 


34.4 


94.7 


106 




0.350 


175 


0.551 


0.244 


33.7 


92.1 


110 




0.400 


210 


0.700 


0.164 


33.3 


90.0 


109 




0.450 


245 


0.852 


0.083 


32.9 


88.2 


111 




0.460 


250 


0.882 


0.066 


32.8 


87.8 


112 




0.470 


255 


0.912 


0.050 


32.7 


87.3 


115 




0.480 


249 


0.941 


0.033 


32.5 


86.5 


110 




0.490 


223 


0.970 


0.015 


35.8 


99.7 


100 




0.495 


210 


0.984 


0.007 


41.7 


117.2 


82 


0.05 


0.202 


55 


0.065 


0.457 


36.4 


102.6 


98 
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Table 5— 


-Continued 








^pvn ( S ) 




S 


X n f 




Z 


A 


n/s 




0.252 


120 


0.236 


0.393 


35.4 


98.3 


104 




0.301 


170 


0.394 


0.320 


34.6 


95.1 


106 




0.351 


217 


0.550 


0.243 


34.0 


92.6 


108 




0.401 


263 


0.702 


0.164 


33.5 


90.4 


110 




0.451 


310 


0.852 


0.082 


33.2 


88.8 


111 




0.461 


315 


0.883 


0.066 


33.0 


88.2 


114 




0.471 


320 


0.912 


0.049 


33.0 


87.9 


113 




0.481 


323 


0.942 


0.033 


32.9 


87.4 


113 




0.490 


295 


0.970 


0.015 


35.4 


97.6 


101 




0.495 


278 


0.984 


0.006 


41.9 


117.7 


83 


0.1 


0.204 


70 


0.067 


0.453 


36.6 


103.2 


97 




0.253 


145 


0.228 


0.386 


35.7 


99.2 


99 




0.302 


210 


0.392 


0.316 


34.9 


95.8 


103 




0.352 


268 


0.548 


0.240 


34.2 


93.2 


106 




0.402 


325 


0.702 


0.162 


33.7 


90.9 


108 




0.451 


385 


0.853 


0.081 


33.4 


89.2 


109 




0.461 


390 


0.883 


0.064 


33.3 


88.8 


108 




471 


400 

t: W \J 


913 


048 


33 2 


88 4 


110 




0.481 


410 


0.944 


0.032 


33.1 


87.9 


113 




0.491 


392 


0.973 


0.015 


32.5 


85.9 


108 




0.496 


369 


0.986 


0.006 


42.5 


118.8 


78 


0.25 


0.210 


105 


0.082 


0.443 


37.1 


104.5 


97 




0.258 


200 


0.239 


0.378 


36.2 


100.2 


99 
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( /~\ttI~ in n on 

vjUII L 111 LlCU. 










S 


X *J 


X n,f 


Z 


A 


n/s 


0.306 


285 


0.401 


0.309 


35.3 


96.8 


103 


0.355 


360 


0.555 


0.234 


34.6 


94.1 


104 


0.404 


435 


0.707 


0.157 


34.1 


91.6 


106 


0.453 


520 


0.859 


0.078 


33.7 


89.7 


110 


0.463 


530 


0.887 


0.061 


33.7 


89.6 


105 


0.473 


544 


0.919 


0.045 


33.5 


88.9 


111 


0.483 


565 


0.949 


0.029 


33.4 


88.4 


115 


0.493 


567 


0.978 


0.012 


33.1 


86.7 


107 


0.498 


532 


0.990 


0.004 


43.7 


121.8 


80 




x/x 




x/x 




x/x 




x/x 




x/x 




x/x 




x/x 




x/x 




x/x 



